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ABSTRACT 

We use new studies of the cosmic evolution of star-forming galaxies to estimate 
the production rate of ionizing photons from hot, massive stars at low and intermedi- 
ate redshifts. The luminosity function of blue galaxies in the Canada-France Redshift 
Survey shows appreciable evolution in the redshift interval z — — 1.3, and generates 
a background intensity at 1 ryd of Jr, m 1.3 x 10 _21 (/ esc ) ergs cm~ 2 s _1 Hz -1 sr _1 
at z w 0.5, where (/ eS c) is the unknown fraction of stellar Lyman-continuum photons 
which can escape into the intergalactic space, and we have assumed that the absorption 
is picket fence-type. We argue that recent upper limits on the Ha surface brightness 
of nearby intergalactic clouds constrain this fraction to be is 20%. The background 
ionizing flux from galaxies can exceed the QSO contribution at z « 0.5 if (/ csc ) 6%. 
We show that, in the general framework of a diffuse background dominated by QSOs 
and/or star- forming galaxies, the cosmological baryon density associated with pho- 
toionized, optically thin gas decreases rapidly with cosmic time. The results of a recent 
Hubble Space Telescope survey of O vi absorption lines in QSO spectra suggest that 
most of this evolution may be due to the bulk heating and collisional ionization of the 
intergalactic medium by supernova events in young galaxy halos. 



1 INTRODUCTION 

The integrated ultraviolet flux arising from quasars and hot, 
massive stars in star-forming galaxies is likely responsible for 
maintaining the high degree of ionization of the intergalactic 
medium (IGM). The large, low-density intergalactic clouds 
which produce the Lyman-a forest lines in the absorption 
spectra of background QSOs represent one of the observa- 
tional signature of this high level of ionization (Bechtold 
1994; Giallongo et al. 1996) . Locally, the ultraviolet ionizing 
background (UVB) may be responsible for the ionization of 
the hydrogen clouds located in the galactic halo (Ferrara 
& Field 1994), and for producing the abrupt truncation in 
the H I distribution at the edge of nearby spiral galaxies 
(Bochkarev & Sunyaev 1977). 

The contribution of QSOs to the UVB is the easiest to 
assess with some degree of reliability. Madau (1992), and, 
more recently, Haardt & Madau (1996) have computed the 
intensity of the UVB as a function of redshift based on the 
most recent estimates of the quasar luminosity function from 
z = to z ~ 5. They included the reprocessing of ultraviolet 
radiation by intergalactic material, and showed that QSO 
absorption- line systems are sources, not just sinks of ionizing 
photons. The resulting metagalactic flux at the Lyman edge 
is found to increase from ~ 10 -23 ergs cm -2 s _1 Hz -1 sr _1 



at the present-epoch to ~ 5xl0 -22 ergs cm -2 s _1 Hz -1 sr _1 
at z — 2.5. The background intensity stays nearly constant 
in the redshift interval 2 = 1.5 — 3.5, to drop rapidly beyond 
z = 4 because of the steep decline of the quasar population. 

Hot, massive stars in star-forming galaxies have also 
been suggested as important contributors to the UVB (Bech- 
told et al. 1987; Songaila, Cowie, & Lilly 1990; Miralda- 
Escude & Ostriker 1990) at early epochs. As the present rate 
of production of metals in normal galaxies is too low to ac- 
count for the observed element abundances, there must have 
been an epoch when the heavy element production rate per 
unit mass was several times larger than it is today (Madau et 
al. 1996). Data on the galaxy population at 2 ^ 3 are sparse 
but are accumulating rapidly (Steidel et al. 1996; Madau et 
al. 1996). The recent detection by Tytler et al. (1995) and 
by Cowie et al. (1995) of C iv in the Lyman-a forest clouds 
has provided the first evidence of widespread chemical en- 
richment in the IGM at z ~ 3. Madau & Shull (1996) have 
computed the ionizing stellar radiation flux which accompa- 
nies the production of metals at high-z, and found that this 
may be significant, comparable to the QSO contribution if a 
fraction 25% of the UV radiation emitted from stars can 
escape into the intergalactic space. At low and intermediate 
redshifts, the Canada- France Redshift Survey (Crampton et 
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al. 1995) has provided new information on the properties 
and evolution of field galaxies at z < 1.3. Finally, the first 
survey for O vi 1032, 1038A absorption lines in QSO spec- 
tra (Buries & Tytler 1996) has shown the likely presence of 
a substantial cosmological mass density of hot, collisionally 
ionized gas at (z) = 0.9. 

In this paper we estimate the contribution of blue, star- 
forming galaxies to the UVB, and use some new upper lim- 
its to the intensity of the UV diffuse radiation field at the 
present-epoch to constrain the average escape fraction of 
Lyman-continuum (Lyc) photons from such systems. We 
also show that the cosmological gas density of T « 20, 000 
K material photoionized by the UVB decreases strongly 
with cosmic time, and suggest that the recent detection of 
O vi intergalactic absorption lines may provide some clues 
on the fate of the "missing" baryons. Throughout this pa- 
per, we shall adopt a flat cosmology with Ho — 50hso km 
s" 1 Mpc" 1 and q = 0.5. 



2 STAR-FORMING GALAXIES AND THE 
UVB AT Z<1.5 




2.1 Blue Galaxy Emissivity 



The results of the Canada-France Redshift Survey have 
made it possible to disantangle the evolution of the luminos- 
ity function (LF) from redshift z — 0.02 to redshift z — 1.3 
of the blue and red galaxy population separately (Lilly et 
al. 1995). At variance with that of red objects, the LF of 
blue galaxies (bluer than present day Sbc) shows significant 
cosmological evolution. This evolution can be represented 
as a brightening of the average luminosity of the blue pop- 
ulation as the redshift increases (although density evolution 
can be present as well). Between z w 0.3 and z w 0.6 the LF 
brightens by about 1 magnitude. Beyond, and up to z « 1.3 
the luminosity evolution of the bright end of the LF levels 
off. For z < 0.2, the new survey shows a significant excess 
relative to the local LF of Loveday et al. (1992). 

In the AB photometric system, the average rest-frame 
ultraviolet color of the blue population, in the redshift inter- 
val considered, is (U — V) ~ 0.8. To estimate the mean num- 
ber of Lyc photons produced by O-B stars, we have modeled 
the intrinsic UV spectral energy distribution of star-forming 
galaxies using the Bruzual & Chariot (1993) stellar popula- 
tion synthesis code. We assume a Salpeter initial mass func- 
tion with lower and upper cutoffs of 0.1 and 125Mq, and 
a constant star formation rate. The fiducial galaxy age is 
fixed at 2 Gyr; this yields the required U — V color af- 
ter convolving with the standard Johnson broadband filters. 
Differences greater than 1 Gyr in the age produce changes 
greater than 0.1 mag in the average colors. 

The comoving emissivity at hvL = 1 ryd of our galaxy 
sample can then be written as 



e(v L ,z) = (/ esc ) 



4>(L, z)LdL, 



(1) 



where (/ CS c) is the fraction of Lyc photons emitted from stars 
which can escape into the intergalactic space, assumed to 
be independent of frequency (picket fence- type absorption), 
in analogy with the leakage of ionizing photons from the 



Figure 1. a) Galaxy emissivity as a function of redshift derived 
from the binned luminosity function of Lilly et al. (1995). b) The 
UV background as a function of redshift produced by galaxies 
(solid thin curve), quasars (dashed curve), and galaxies + quasars 
(solid thick curve). The galaxy contribution has been computed 
assuming an escape fraction of ionizing photons into the IGM of 
(Use) = 15%. 



galactic disk. E v is the spectral energy distribution. The 
luminosity function in the B-band, 4>, at various z is taken 
from Lilly et al. (1995), and we have adopted a constant 
value of Lmin, corresponding to Mb = —17.8. The ensuing 
galaxy emissivity at the Lyman edge is plotted in Figure 1 
as a function of redshift: a strong cosmological evolution is 
apparent from z = to z fa 1. 

Our proce dure yields e(912A, 0) = 2 x 10 25 {/ csc ) ergs 
s~ Hz - Mpc . To check for systematic errors, we have re- 
computed the present-day ionizing emissivity starting from 
the Ha luminosity density of the local universe, 1.3 ± 0.7 x 



10' 



59 ergs s 1 Mpc 



(Gallego et al. 1995). Assuming case- 
B recombination theory, an escape fraction of 50%, and our 
fiducial population synthesis galaxy spectrum, this value im- 
plies e(912A, 0) = 1.1 ± 0.5 x 10 25 ergs s" 1 Hz" 1 Mpc" 3 , in 
good agremeent, within the errors, with the extrapolated 
value plotted in Figure 1. 



2.2 Ionizing Background Intensity 

The mean specific intensity J„ of the UV background at the 
observed frequency u , as seen by an observer at redshift z , 
is given by 

(l + z ) 3 



J{v ,z ) 



-f 



dz 



dz (1 + zf 



e(u, z) exp[— r e //(i/ , z , z)\(2) 



where r e // is the effective photoelectric optical depth of a 
cloudy intergalactic medium following Madau (1992), e(v, z) 
is the proper volume emissivity at frequency v — v a (\ + 
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z)/(l + z ), and dl/dz is the line element in a Friedmann 
cosmology. 

The main uncertainty in the estimate of the UVB gen- 
erated by star-forming galaxies is the escape fraction of Lyc 
photons (/esc) into the IGM. While theoretical models of the 
structure of H n regions in a vertically stratified disk suggest 
an escape fraction of ~ 15% through the H i layers into the 
halo (Dove & Shull 1994), the recent data obtained by the 
Hopkins Ultraviolet Telescope show that less than 3% of the 
stellar UV radiation can escape from the four, low-redshift 
starburst galaxies observed by Leitherer et al. (1995). Given 
the peculiar nature of the selected targets, however, it is 
conceivable that the escape fraction might be significantly 
reduced in these objects by the presence of dust absorption. 
From a reevaluation of the same data Hurwitz et al. (1997) 
suggest less restrictive upper limits with an average value 
closer to 19%. 

An indirect argument, put forward by Patel & Wilson 
(1995a, b), suggests a Lyc escape fraction from the galactic 
disk into the halo greater than 50%. From their observations 
of Ha emission in nearby galaxies, these authors were able 
to estimate the total star formation rate and compare it 
with the value derived from the distribution of OB stars 
in the galactic disk. In the presence of ionization balance, 
the two values should be consistent. On the contrary, heavy 
losses of ionizing photons appear to be taking place in the 
field. In this respect, it is important to note that, although 
the fraction of ionizing photons leaking from the galactic 
disk could be relatively high, Lyc absorption by cold gas 
in galaxy halos may be quite common, as suggested by the 
identification of the galaxies associated with the Lyman- 
limit systems observed in quasar spectra (e.g., Steidel 1995). 

The study of nearby isolated intergalactic clouds where 
no sign of stellar activity is detected can provide a stringent 
upper limit to the local UVB and hence to the mean es- 
cape fraction. In particular, by looking at the Ha radiation 
that should escape from these optically thick clouds due to 
the reprocessing of the incident Lyc photons, Vogel et al. 
(1995) (see also Donahue, Aldering & Stocke 1995) have re- 
cently set a 3-ct upper limit of Jl < 8 x 1CT 23 ergs cm" 2 s" 1 
Hz -1 sr" 1 at z = 0. Given the constraints on the luminosity 
function and the average spectral shape of the blue galaxies, 
the local UVB limit implies that no more than 20% of Lyc 
photons can escape into the IGM from star-forming regions. 
The evolution of a galaxy-dominated UVB is shown as a 
function of redshift in Figure 1. Adding together the galaxy 
contribution with the one estimated from quasars and lo- 
cal AGNs (Haardt & Madau 1996), we find that a value of 
(/esc) = 15%, predicted by theoretical models (Dove & Shull 
1994), yields a total UVB which is still consistent with the 
local upper limits. It is interesting to notice that, since the 
opacity of the IGM at z < 1.5 is low, galaxies at intermedi- 
ate redshifts can provide a non negligible contribution to the 
local metagalactic flux. The current limits on the UVB from 
Ha-brightness studies in fact provide the most stringent con- 
straints to the total ionizing emissivity at z ~ 0.5 — 1, where 
the bulk of the blue galaxy evolution is observed. 

As shown in Figure 1, an escape fraction of (/esc) = 15% 
implies a total UVB of Jl ~ 10~ 22 ergs cm" 2 s" 1 Hz -1 
sr" 1 at 2 ~ 0.15. If the contribution of the blue galaxies to 
the metagalactic flux extends to z ^ 1.5, i.e. if the escape 
fraction of Lyc photons into the IGM is significant even at 



early epochs, it is possible to envisage a scenario where the 
background intensity remains nearly constant, Jl = 2 — 4 x 
10" 22 ergs cm" 2 s" 1 Hz -1 sr" 1 , from z = 0.4 to z = 4, i.e. 
for a fraction of 50% of the age of the universe. Alternatively, 
QSOs will dominate the metagalactic flux if (/ eS c) < 6%. 
The intensity, spectrum, and evolutionary history of a QSO- 
dominated UVB have been recently computed by Haardt & 
Madau (1996) on the basis of our current knowledge of the 
QSO luminosity function. The integrated diffuse flux at the 
Lyman edge is found to increase from 2 x 10" 23 at z — to 
k5x 10" 22 ergs cm" 2 s" 1 Hz -1 sr" 1 at z « 2. Over the 
redshift range z = 2 — 3.5, Jl remains roughly constant, to 
decrease to 2 x 10" 23 ergs cm" 2 s" 1 Hz -1 sr" 1 at z w 5. 

Since the luminosity function of star forming galaxies at 
z > 1 is poorly known it is difficult to estimate the galaxy 
contribution to the UVB in the redshift interval 2 < z < 5. 
However an extrapolation of the z ~ 1 luminosity function 
at higher z would imply a total UVB at z ~ 3.5 of the order 
of J L =6±1 x 10" 22 ergs cm" 2 s" 1 Hz" 1 sr" 1 . 

An independent estimate of Jl at high redshifts comes 
from the statistical study of the proximity effect in the spec- 
tra of QSOs (Bajtlik, Duncan, & Ostriker 1988; Bechtold 
1994). The value recently derived from a large high resolu- 
tion sample of Lyman-a forest lines is Jl — 5 ± 1 x 10" 22 
ergs cm" 2 s" 1 Hz -1 sr" 1 , constant in the redshift range 
1.7 < z < 4.1 (Giallongo et al. 1996). As a direct conse- 
quence, while in a QSO-dominated background model the 
photoionization rate will remain approximately constant 
from 2 ~ 3.5 to z ~ 1.5, to drop by a factor of ~ 30 by 
the present epoch, in a universe in which star-forming galax- 
ies contribute significantly to the metagalactic flux, Jl may 
only decrease by a factor ~ 6 from z — 2 to the present 
epoch. 



3 THE BARYON DENSITY OF 
PHOTOIONIZED GAS 

A simple estimate of the cosmological density of the baryons 
hidden in the Lyman-a forest clouds can be obtained from 
the following quantities: the intensity Jl of the UVB, the 
number density of lines along the line of sight dN/dz, the 
typical cloud radius R, and the observed H i column. The 
latter spans a large range of values, from 10 12 cm" 2 to 10 15 
cm" 2 and more, showing a bending below Nm « 10 14 cm" 2 
where the slope of the power-law distribution flattens from 
1.8 to 1.4 at (z) = 3 (Giallongo et al. 1996). The density of 
lines with Nm ^ 10 14 cm" 2 , as derived from Hubble Space 
Telescope (Bahcall et al. 1996) and optical data (Giallongo 
et al. 1996), stays approximately constant - dN/dz ~ 25 - 
from z ~ 0.15 to z ~ 1.8, to increase rapidly for z > 2, 
dN/dz ~ 220 at z « 3.8. 

The baryon density of photoionized gas can be ex- 
pressed in units of the critical value Higm = 
fvriH/n CT it where fv is the volume filling factor and riH 
is the typical hydrogen gas density of individual clouds. 
The volume filling factor is proportional to the number of 
systems along the line of sight and depends on the cloud 
geometry, while the gas density is a measure of the pho- 
toionization state of the absorbers (Madau & Shull 1996), 
n H ex R- 1/2 N^ i 2 T a363 J 1/2 (a cos(9) 1/2 . The temperature 
dependence is very weak for highly photoionized clouds: a 
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value of T — 2 x 10 4 K has been adopted here. The as- 
pect ratio a = R/l generalizes the absorber geometry from 
spheres to disks of transverse radii R to half-thickness I. For 
(a cos 6) = 1 the usual limit of spherical clouds is obtained. 

Hence, the cosmological mass density produced by the 
Lya clouds is proportional to 

O — 1 q v in-H ta+3\- 1 / 2 ( J-22-R100 \ ^pO.363 (r>\ 

52igm - 1.3 x 10 [—) ^ - cosfl j T4.3 (6) 

where J = J_22l0~ 22 ergs s _1 cm -2 Hz -1 sr _1 , R = 
flioolOO kpc, T = T 4 . 3 2 x 10 4 K, a = 1.5 is the spectral 
slope of the ionizing UVB and B(z) is the number of lines 
per unit Nhi and unit redshift interval derived at given 
redshifts from the observed spectra as in Giallongo et al. 
(1996). Thus, dN/dz = B(z) J N~fdN m . We include the 
contribution from Lya lines of various column densities on 
the basis of the known shape of the column density distribu- 
tion. For JVhi > 10 14 cm~ 2 this distribution is fairly steep, 
with /3 S ~ 1.8, and the contribution to JIigm becomes pro- 
gressively small. For Nm < 10 14 cm -2 , [3f ~ 1.4 and the 
contribution to the mass density parameter increases slowly 
with decreasing JVhi. We have adopted a lower Nm cutoff 
of 10 12 cm -2 as derived from high resolution data (Webb et 
al. 1992, Giallongo et al. 1995, Hu et al. 1995). 

The cloud size and geometry are subject to substantial 
uncertainties. An estimate of the characteristic size of ab- 
sorbers is provided by the statistical coincidence of absorp- 
tion lines in closely separated quasar pairs. Recent obser- 
vations of the quasar pair 1343+264A/B at z ~ 1.8 (Bech- 
told et al. 1994, Dinshaw et al. 1994) have shown that the 
Lya sizes are of the order of R ~ 200h 50 1 kpc, much larger 
than previously thought. Observations at lower redshifts, 
z ~ 0.5, show even larger sizes, R ^ 360h | T 1 kpc, indepen- 
dently of the cloud structure or geometry (Dinshaw et al. 
1995). The mild redshift evolution implied by these prelim- 
inary measurements is consistent with the expectations of 
the standard CDM cosmological scenarios. As outlined by 
the numerical simulations of Miralda-Escude et al. (1996), 
although the gas moves on average to regions of higher over- 
density, the dominant effect for the evolution of the average 
cloud properties appears to be the Hubble expansion. 

For these reasons, we have assumed a typical radius 
R — 200h 50 1 kpc at z = 1.8 which evolves in time following 
the Hubble expansion as R oc (1 + z) -1 . 

The ensuing cosmological baryon density due to Lyman- 
q clouds is shown in Figure 2 from z = to z ~ 4 as 
a function of cosmic time. The redshift dependence of the 
baryon density of photoionized gas is associated with the 
product {JR) l t 2 dN /dz where the dominant factor is the 
number density evolution of the clouds. 

Two different evolutionary scenarios can be envisaged 
at this point, depending on the sources of the ionizing back- 
ground. In Figure 2a Qigm(^) has been computed for a 
QSO-dominated UVB. In Figure 2b both galaxies and QSOs 
contribute to the metagalactic flux. In the redshift inter- 
val 1.7 < z < 4.1, the UVB assumes the constant value 
J912 ~ 5 x 10~ 22 ergs cm~ 2 s~ x Hz -1 sr _1 , as derived 
from the proximity effect. The increase of the baryon den- 
sity at higher redshifts is due to the increase of the num- 
ber density of Lyman-a lines, which is partly compensated 
by the decrease of the average cloud sizes. Lyman-a ab- 
sorbers at z ?s 4 can easily account for all the baryons in 
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Figure 2. Mass density parameter of photoionized intergalactic 
gas as a function of cosmic time, a) QSO-dominated UVB. b) 
Same including a contribution from star-forming galaxies. Points 
are derived from the Lyman-a sample described in the text. Two 
exponential curves are shown simply for comparison (see sect. 4), 
with e-folding timcscalcs equal to 1 and 8 Gyr. Again, the galaxy 
contribution has been computed assuming an escape fraction of 
ionizing photons into the IGM of (/esc) = 15%. 

the universe predicted by nucleosynthesis, f^lo = 0.05 ± 1 
(Walker et al. 1991). Note that a large aspect ratio, a ^ 10, 
must be adopted to avoid Oigm values significantly in excess 
of the nucleosynthesis constraints, as outlined by Rauch & 
Haehnelt (1995). 

The residual baryon density is already about 30% at 
z ~ 1.8 (i.e., t ~ 3 Gyr). By contrast, for z < 1.7 the number 
density of Lyman-a clouds stays nearly constant and the 
evolution of the density parameter depends mainly on the 
evolution of the ionizing sources contributing to the UVB 
and of the cloud size. In any case, a further decrease of Oigm 
is present in the redshift interval z = 0.3 — 1.7. This evolution 
is stronger in the case of a QSO-dominated UVB where a 
value of f2iGM — 0.007 is found at z = 0.3. In a galaxy- 
dominated UVB the evolution is slower and Oigm — 0.01 at 
the same redshift. These values are larger than the baryonic 
"visible" mass density of galaxies, f2» = 0.004e ±0 ' 3 (e.g., 
Peebles 1993). 



4 DISCUSSION 

We have presented here two main results. First, we have 
computed the integrated emission of ionizing photons from 
the observed population of star-forming galaxies at z < 1.3, 
and shown that it may equal or even exceed the overall QSO 
contribution to the UVB in this redshift range if the mean 
fraction of UV photons which can escape from individual 
galaxies into the IGM is not negligible. We have shown in 
§ 2 that (/esc) is constrained to be < 20%, in order to satisfy 
the upper limit on the intensity of the local UVB. If we adopt 
(/esc) ~ 15% as a fiducial value for the escape fraction, the 
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ensuing UVB is found to evolve little in the redshift interval 
between 2 = 0.4 and z = 4. 

We have also shown that the baryon fraction which 
is associated with the photoionized Lya clouds appears to 
decrease rapidly with cosmic time. At 2 ~ 4, Hiqm may 
account for all of the nucleosynthesis baryons of the uni- 
verse. As shown in § 3, the mass density parameter de- 
pends only weakly on the intensity, Jl, of the UVB, with 
fiiGM oc J 1 / 2 . Hence tighter constraints on the escape frac- 
tion from galaxies are unlikely to change qualitatively the 
evolution of Oigm (see Fig. 2). The dependence of J7igm on 
the typical cloud radius and geometry is also relatively weak, 
f2iGM oc (R/a) 1 / 2 = I 1 / 2 . Moreover, in CDM scenarios, lines 
with iVm < 10 15 cm -2 have typical sizes which appear es- 
sentially uncorrelated with the H i column (Miralda-Escude 
et al. 1996). Thus, while the exact normalization of Higm 
depends on these values, it is difficult to explain the strong 
evolution inferred in the redshift interval 2 = 1. 7 — 4 as a 
simple geometrical effect. Note that photoionization mod- 
els and the statistics of coincident absorptions in QSO pairs 
suggest a typical cloud thickness in the range 21 « 100 — 300 
kpc, independently of the assumed aspect ratio (Fang et al. 
1996). At z ~ 4, the baryon content of the Lyman-a forest 
clouds may then be significant both in the case of spherical 
or disk-like geometry. 

We can parameterize the cosmological evolution of 
fiiGM as an exponentially decreasing function with cosmic 
time, f2iGM oc exp(— t/r). As shown in Figure 2, this yields 
two characteristic e-folding timescales, r ~ 1 Gyr for t < 3 
Gyr and r ~ 8 Gyr at later epochs. Thus, roughly 60% of 
the photoionized gas in the universe "disappears" at early 
epochs over a timescale which is much shorter that the Hub- 
ble time. What is the fate of these "missing" baryons? They 
are unlikely to fragment into stars, as the baryonic visible 
mass density of galaxies today is known to be only a small 
fraction of the value derived from cosmological nucleosynthe- 
sis. The decrease of SIigm observed in the Lya line popula- 
tion must mainly reflects either a dilution of the photoion- 
ized gas clouds in a general diffuse intergalactic medium or 
an increase of the fraction of the IGM which is collisionally 
ionized at temperatures T ^ 10 5 K by supernovae shocks. 
If a strong evolution in the neutral column density of clouds 
were present following the Hubble flow, it would produce 
at 2 = 1.7 lines ten times weaker {[ogNni ~ 11) than at 
2 ~ 3.8. However, since the contribution to Qigm increases 
slowly for decreasing Nhi {Qigm oc Nhi for /3 ~ 1.4), lines 
ten times weaker in Nhi increase Q. only by ~ 0.005 when 
2 decreases from z ~ 3.8 to z ~ 1.7 (i.e. from t ~ 1 Gyr 
to t ~ 3 Gyr). Thus the strong f2 evolution derived in the 
redshift interval z = 1.7 — 4 can not be mainly due to a 
decrease (with decreasing redshift) of the lower cutoff of the 
Nhi distribution unless the sizes evolved much faster than 
assumed in this simple model. 

A plausible possibility appears to be the additional 
heating of the absorbing gas at temperatures T ~ 10 5-6 K 
by the gravitational accretion into progressively more mas- 
sive halos, with higher velocity dispersions, or by collisional 
ionization from supernovae winds. The easiest way for find- 
ing collisionally ionized, cosmologically distributed material 
at T ^ 10 5 ' 5 K is to search for O vi absorption. O vi is 
most prevalent at these temperatures, while O vi lines are 
stronger than those of C iv for T ^ 10 6 K. The recent results 



of the first survey for O vi 1032, 1038A absorption lines in 
QSO spectra (Buries & Tytler 1996) suggest the presence of 
a substantial cosmological mass density of hot, collisionally 
ionized gas at (z) = 0.9. If the bulk heating were mainly due 
to supernovae explosions in spheroidal systems, as suggested 
by recent numerical simulations (e.g., Miralda-Escude et al. 
1996), the strong evolution of Oigm observed between 2 = 2 
and 2 = 4 could be triggered by the star-formation activity 
in galaxies at high redshift. Note that, while an e-folding 
timescale of only 1 Gyr is much shorter than the decay time 
of star formation, w 4 — 7 Gyr, characteristic of late type 
spirals, it is comparable with the decay time of star forma- 
tion characteristic of "average" quiescent elliptical galaxies 
(e.g., Bruzual & Chariot 1993). Various authors have argued 
in favour of a significant contribution to the UVB by star- 
forming galaxies (Bechtold et al. 1987; Miralda-Escude & 
Ostriker 1990; Madau 1991; Madau & Shull 1996). Order-of- 
magnitude arguments suggest that the overall kinetic energy 
released in supernova explosions can be about one-third of 
the radiative UV power (Miralda-Escude & Ostriker 1990). 
In this case a highly photoionized IGM at T w 10 4 ' 5 K may 
be quickly heated to temperatures 10 times larger (Giroux 
& Shapiro 1996). If this scenario will be confirmed by new 
observational constraints on R and Jl and by detailed the- 
oretical models, then the evolution of the Lyman-a forest 
clouds may be used as a probe of the cosmic star formation 
rate as a function of time. 
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